The Phonon Drag Effect in Single- Walled Carbon Nanotubes 
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A variational solution of the coupled electron-phonon 
Boltzmann equations is used to calculate the phonon drag 
contribution to the thermopower in a 1-D system. A sim- 
ple formula is derived for the temperature dependence of 
the phonon drag in metallic, single- walled carbon nanotubes. 
Scattering between different electronic bands yields nonzero 
values for the phonon drag as the Fermi level varies. 
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I. INTRODUCTION 

Since their discovery [1], carbon nanotubes have pro- 
vided a testbed for fundamental and applied physics. 
Characterization of these new systems is crucial for future 
progress. The thermopower of a material is an intrinsic 
quantity which yields important information related to 
the electronic band structure, electron-phonon coupling 
parameters, and relaxation rates of the system. Recent 
experiments on mats of single walled carbon nanotubes 
[2-8] have found surprisingly large values for the thermo- 
electric power (~ 50/jV/K) under various ambient condi- 
tions. When the thermopower is plotted versus temper- 
ature most of the data show a quasi-linear behavior at 
temperatures greater than 200K. Thermopower linear in 
temperature suggests that conduction through metallic 
tubes dominates the thermopower at these temperatures. 
However, some data show nonlinearities below ~ 200K . 
Here a peak in the thermopower is observed. The origin 
of this peak is the subject of much debate. Several effects 
that could lead to such behavior have been discussed in 
the literature. These include parallel transport through 
semiconducting tubes [2] , a one-dimensional Kondo effect 
[4], and the phonon drag effect. Semiconducting tubes 
are not expected to contribute much to the thermopower 
because metallic tubes have a larger electrical conductiv- 
ity [3]. The Kondo effect has been invoked to explain 
the large peak in samples containing magnetic impuri- 
ties. However, a smaller, broad peak remains when the 
magnetic impurities are removed and therefore does not 
explain the universality of the low temperature nonlin- 
earity in samples without magnetic impurities. 

Turning to the phonon drag effect, very little theoret- 
ical work has been done to calculate the phonon drag 
contribution to the thermopower in metallic carbon nan- 
otubes. Rough, low temperature, estimates often quote 
the result 



S, 



drag 



OC C v 



(1) 



where C v is the lattice specific heat. However, this for- 
mula assumes that phonons scatter electrons within one 
parabolic band and that electron-phonon scattering acts 
as the dominant phonon decay mechanism. Applying 
Eq. (1) to nanotubes with two bands, particle and hole, 
yields a negligible contribution to the thermopower due 
to drag because the particle and hole contributions can- 
cel. 

The objective of this work is to extract the tempera- 
ture dependence of the phonon drag contribution to the 
thermopower for a one dimensional system, within the 
linear band approximation. We consider a model wherein 
mechanisms other than electron-phonon scattering limit 
phonon lifetimes. We derive a new expression for the 
phonon drag in these systems using a solution of the cou- 
pled electron-phonon Boltzmann equations. We find that 
mterband scattering gives a non-zero contribution to the 
thermopower when the Fermi level does not lie at the 
band crossing. Interband transitions near the Fermi level 
contribute appreciably to the phonon drag thermopower. 
Our results can be summarized in the following formula 



Stotal = AT - 



in the limit 
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where A and B are fitting parameters, \i is the energy 
difference between the Fermi level and the band cross- 
ing, v the electron speed, c a typical phonon speed, and 
To the Debye temperature, m parameterizes the tem- 
perature dependence of the phonon relaxation time. In 
metals one typically finds m ~ 1 at temperatures near the 
Debye temperature. The term linear in temperature is 
the usual diffusive contribution for metals and the second 
term is our result for the phonon drag part at low temper- 
atures. This formula yields a peak, due to phonon drag. 
We therefore associate the phonon drag effect with the 
low temperature nonlinearities observed in thermopower 
measurements on single walled carbon nanotubes. 

The plan of the paper is as follows. In Section II we 
outline the basic theory of thermopower measurements 
in metals. In Section III Bailyn's theory of phonon drag 
in metals is reviewed. In Section IV we apply Bailyn's 
theory to a simple 1-D model of single walled, (10,10) 
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armchair carbon nanotubes. Here we derive an expres- 
sion for the phonon drag in these systems. In Section V 
we extract a simple formula by making a low temperature 
approximation. In Section VI we conclude by discussing 
the limits of our approximations. 



II. THERMOPOWER DUE TO ELECTRON 
DIFFUSION AND PHONON DRAG 

The application of a temperature gradient to a metal 
leads to the diffusion of charge carriers from the warm 
to the cold end of the sample. The thermopower mea- 
sures the charge build up across the sample. It is given 
by S = lim A T^o AV/AT, where AV and AT arc the 
potential and temperature differences across the sample, 
respectively. The thermopower due to diffusion, Sdiff, can 
be calculated from the standard set of transport coeffi- 
cients. The transport coefficients can in turn be derived 
from the Boltzmann equation for the electron distribu- 
tion function. A general argument leads to the Mott 
expression for the thermopower due to diffusion 



k%T ( diner \ 



3 |e| \ dE J 



(3) 



EF 



where E is the energy, a the electrical conductivity, and 
EF the Fermi energy. When there is more than one band 
present the thermopower due to each band adds 



Sdifi 



J2e g^dig 



(4) 



where i is the band index. 

We apply the above equations to the case of metallic 
carbon nanotubes by considering electrons in one dimen- 
sion. The electronic states fill two overlapping, parabolic 
bands, the particle and hole bands, up to the Fermi level. 
With only one parabolic band we expect the diffusive 
part of the thermopower in a metal to vary linearly with 
temperature. The factor (^pf) EF depends on the de- 
tails of the system, including the density of states. How- 
ever, with two bands present, the thermopower due to 
states filling the hole band cancels the contribution from 
the particle band when the Fermi level lies at the band 
crossing. Within this approximation we have no net con- 
tribution to the thermopower in a metallic tube. 

When the Fermi level is allowed to move within the 
rigid 7T bands the resulting thermopower is non-zero [3]. 
Recent calculations show that an enhancement in the 
density of states due to impurities [3,9,10] or tube-tube 
interactions [2,11,12] may generate large contributions to 
the thermopower through the term {^jf) EF - We there- 
fore consider the following standard form for the diffusive 
contribution to the thermopower in metals 



where the constant A is a fitting parameter which may 
vary with the Fermi level. Recent experiments on mats of 
single walled, carbon nanotube bundles do indeed show 
this behavior at large temperatures (T > 200 A") [3]. 

An anomalous peak in the thermopower appears at low 
temperatures in several different experiments on single- 
walled and multi-walled carbon nanotubes. In standard 
thermopower measurements of 3-D metals such behavior 
is often associated with the phonon drag effect whereby 
the phonon flux from the hot end of the sample to the 
cold end drags additional charge carriers to the cold end 
of the sample via momentum transfer. This effect adds 
to the thermopower in conventional metals 



Stotal — Sdiff + Sd 



rag 
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Standard, low temperature estimates of the phonon drag 
contribution to the thermopower rely on the relation 



^drag 
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(7) 
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where C v is the lattice specific heat, n the carrier density, 
and the t is the transfer factor, t is a rough estimate of the 
probability that a phonon collides with an electron rela- 
tive to all scattering events. In one dimension the lattice 
specific heat is nearly linear in temperature [14,15] and 
does not provide the nonlinear temperature dependence 
required to explain the peak observed in measurements 
on carbon nanotubes [16]. Moreover, if we consider only 
intraband scattering, the contributions from states filling 
the electron and hole bands should cancel to give no net 
drag. 

The above formula for the drag contribution cannot 
be applied to metallic carbon nanotubes for two rea- 
sons. The derivation of the above formula [18] relics 
on a free, electronic band structure where transitions lie 
only within the parabolic bands. It also assumes that 
the dominant decay mechanism for phonons is electron- 
phonon scattering. Below we derive a new formula for the 
low temperature phonon drag contribution to the ther- 
mopower in one dimension. We find a nonzero contribu- 
tion to the phonon drag part of the thermopower when 
we include transitions between two linear bands and as- 
sume that the dominant decay mechanism for phonons is 
not electron-phonon scattering. 



III. BAILYN FORMALISM FOR PHONON DRAG 

The phonon drag contribution to the thermopower 
in a metal may be calculated by solving the coupled 
electron-phonon Boltzmann equations. In this section 
we briefly review Bailyn's formalism for calculating the 
phonon drag. Following Ref. [19] we write the Boltzmann 
equation for the electron distribution function, /, in the 
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relaxation time approximation. Using first oder pertur- 
bation theory for the electron transition probabilities one 
finds 

df 



( "^~) co11 ~ E C k,k>j{ S (-) S k>,k+q 



k't' 



[ -N(qj)f(l - /') + (N(qj) + 1)/'(1 - /)] 
+ N(qj)f'(l-f)}} 

-(/-/o)/r(fc) (8) 

where N(qj) is the phonon distribution function, = c\q\ 
is the frequency of a phonon with speed c, /o is the Fermi- 
Dirac distribution function, and r(fc) is the relaxation 
time due to the electron-electron interaction. We define 
the energy delta functions 



S(±) = 5{E{k'i') - E{k£) ± tvjj) 



(9) 



The factor g, . is related to the electron-phonon matrix 
elements 



Mm 



k < k 'i 2fru(qj) 
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where M(qj) = \ < k'£'\VU ■ e{qj)\k£ > \ 2 /MN c is the 
square of the matrix element for the scattering of an elec- 
tron from wave vector k and band £ to wave vector k' and 
band £' by a phonon of wave vector q and polarization j. 
Here VZ7 is the gradient of the ion potential, N c the num- 
ber of cells in the periodic block, M the ion Mass, and 
e(qj) the phonon polarization vector. The above matrix 
element ignores Umklapp scattering [20] . 

From the above expression for the electron distribu- 
tion function we can read off the necessary terms for 
the phonon Boltzmann equation. The second and third 
terms show gains in the phonon distribution. The first 
and fourth terms show losses in the phonon distribution 
function. This gives 

cwm _ ( dNm 

dt -' drift ~ ^ qi > cca 

k,k> 

+ (N(qj) + 1)/'(1 - /)] - 5{+)5 % , rk+ . 
[ -(N(qj) + 1)/(1 - /') + N(qj)f(l /)]} 
- (N(qj) - N )/r(q) (11) 

where r(q) is the phonon relaxation time and N$ is the 
Bose distribution. The above two equations for the elec- 
tron and phonon distribution functions can be solved us- 
ing a variational procedure. From the relevant transport 
coefficients the phonon drag contribution to the ther- 
mopower can be extracted. The most general form for 
which was derived in Ref. [21]. 



^drag 
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x J2 ^r^ l ^ l ')^U T u"^i' T w^-%3 (12) 
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Here a is the electrical conductivity, d the dimensionality 
of the system, the 2 results from a sum over the spin 
degrees of freedom, is the electron group velocity, 
and Vjfj is the phonon group velocity. The factor a is 
the relative probability that the qj phonon will scatter 
an electron from the state k£ to the state k'f, relative to 
all other possible phonon collisions. 

Details of the electron-phonon interaction are included 
in a. Symbolically 



ep 



E 



+ T p 



(13) 



where r ep is the phonon relaxation time due to the 
electron-phonon interaction and r p is the phonon relax- 
ation time due to any other interaction. This may include 
phonon-phonon, phonon-boundary, phonon-impurity, or 
phonon-defect scattering. Using the above results for the 
phonon relaxation rates Bailyn finds [22] 



a{qj;k£,k'£')- - T 9Nq 



^kl,k't' 



T Cll\Q I r 

hLUT p (q) dT "T l^kl,k't ± kl,k't 



(14) 



where 



7 ^"' = mJk^f fa{E{k)){1 /o(i?( ^' ))} 

No(qj)M(qj)5(-)8p U? (15) 

Equations (12), (14), and (15) constitute Bailyn's theory 
of phonon drag in metals. The remainder of this article 
will be concerned with the application of this formalism 
to the case of metallic carbon nanotubes. 



IV. BAILYN FORMALISM APPLIED TO 
METALLIC CARBON NANOTUBES 

We consider a one dimensional lattice lying on the z 
axis with left and right moving electrons with the band 
structure of a (10,10) armchair carbon nanotube. It is 
assumed that the electron relaxation time has a weak 
wave vector dependence for transitions about the K point 
(27r/3a for lattice spacing a) so that r(k£ + 27r/3a) ~ 
T(27r/3a). Our expression for the phonon drag part of 
the thermopower then reads 



S ( i 



rag 



2|e|fc B r v dN (gj) 

QJ 



E 

kl-k'l' 



a(q,j;k£,k'£')[vj; e -v Pe ,}-V qj (16) 
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We note in passing that ignoring all but electron-phonon 
scattering reduces the drag to Eq. (1) when we take a = 
5p [18]. The resulting sum over k, k' becomes hw. We 
assert that this limit is not applicable to metallic carbon 
nanotubes. 

To make progress we evaluate Eq. (16) using approx- 
imations valid for metallic carbon nanotubes. In these 
systems it is reasonable to assume that at small wave 
vectors and energies only acoustic phonons scatter elec- 
trons and that their dispersion is linear [23]. We then 
have 



V qj = cjsign(q)z 



(17) 



where Cj is the phonon speed in the jth band. The elec- 
tronic band structure consists of two nearly linear bands 
crossing at zero energy, Fig. 1. This approximation holds 
as long as the Fermi level lies no more than ~ leV from 
the band crossing. Above leV other bands will compli- 
cate the spectrum. For carrier excitations around the K 
point we have 



E(k£) = ihvk 



(18) 



where v is the Fermi speed and I = +1,-1 labels the two 
bands. With the linear band approximation the electron 
group velocity is then 



(19) 



Substituting the electron and phonon velocities into the 
expression for Sdrag gives 



2\e\k B rv ^ . dN (qj) 
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dT 



(20) 



ki-k' 



Note that the above expression for the phonon drag 
vanishes when only intrah&nd scattering is allowed, i.e. 
1 = 1'. In what follows we consider only znierband scat- 
tering, I' . Define the transfer factor to be 



t(q) = J2 a(qj;k£,k't)[e-t] 

U;k'l' 



(21) 



where a is given by Eqs. (14) and (15). The drag then 
has the simple form 



2\e\kBTV 



^c 3 sign(g) 



dN Q (qj) 
dT 



t(q) (22) 



To simplify the expression for t(q) we make an assump- 
tion regarding the available phonon scattering processes. 
Note that phonons traveling along a single nanotube may 
scatter through many different mechanisms. The large 
amount of surface area exposes the phonons to impuri- 
ties, defects, and neighboring tubes. Furthermore, small 
tube lengths, ~ 10/ito, allow for phonon-boundary scat- 
tering. We consider here two simple forms for the phonon 
relaxation time, valid at low and high temperatures, 
respectively. At low temperatures, the phonon decay 
mechanism is a competition between boundary and de- 
fect scattering which add no temperature dependence to 
the phonon relaxation time. At high temperatures, near 
room temperature, one expects phonon-phonon scatter- 
ing to contribute significantly to phonon decay. Anhar- 
monic scattering will require three phonons. The third 
being an optical mode at long wavelengths. The scatter- 
ing rate will be proportional to the number of available 
optical phonons. One can then show that the relaxation 
time is inversely proportional to temperature, at large 
temperatures. We therefore take the following form for 
the phonon relaxation time 



(24) 



T o [ if 



where ^ is a characteristic temperature, m = for bound- 
ary/defect scattering, and m = 1 for phonon-phonon 
scattering. Comparison with experiment will require a 
more accurate expression for t p . 

We now invoke the assumption that mechanisms other 
than electron-phonon scattering limit the phonon life- 
time, i.e. t p is small [24]. More precisely 



T dNp 
fuvT p (q) ~df~ 



ki.k'f 



(25) 



The above approximation enters phonon drag studies 
of quantum wires in GaAs. Analogous results for the 
phonon drag were also obtained in Ref. [25]. The trans- 
fer factor becomes 



t(q) 



2hjTp(q) 

rr dNp 
1 dT 

E 

fe,l;fc',- 



h. 



Uk' 



y ^ h,-l;k',l 
fe,-l;fe',l 



(26) 



To further evaluate the transfer factor we return to the 
matrix element in Eq. (15). Working with the deformablc 
ion model at low wave vectors (q <C ^Dcbyc) one finds [21] 



To find t(q) we need to evaluate the following expres- 
sion 

E «[* - *] = 2 — fc ; /M ^r: ^ kk ' 7 " — ( 2 3) 

kl,k'l' huT p (q) dT + t-ikiS't 1 1 ki,k'i' 



M(qj) = D^hjjj 



(27) 



where Dj is a constant depending on the ion mass, the 
deformation energy, and other tube parameters including 
the radius and lattice spacing. This constant has been 
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evaluated for (10,10) carbon nanotubes [26]. The low 
wave vector approximation made above has been moti- 
vated by measurements of the Debye temperature. It has 
been measured to be near 1000K in metallic carbon nan- 
otubes [17]. It will be shown that, at low temperatures, 
the temperature dependence of the drag does not depend 
on the precise form of the matrix elements M.. 

The transfer factor contains two terms of the form 

v j = M(qj) 
2^k,k< ki,k>t> 2k B Thuj 

x ^f (E(k£)){l - f (E(k'£'))}N (qj) 

k,k' 

X S(—)5k',k+q 

= M(qj) 
2k B Thuj 

x J2f (E(ke)){l-f (E(k + q,£'))}N (qj) 

k 

x S(E(k + q, £') - E(k£) - hw) (28) 

Where the Fermi-Dirac distribution function is fo{E) = 
j e (B-^)/fe B T + 1 j-i_ The energy delta function in Eq. (28) 
needs to be evaluated for two cases. The first case is £ = 1 
and £' = —1. We then have 



E(k + q,-l)-E(k,l)-fko 
= h[-v(k + q)-vk- c\q\] = 



(29) 



This has a solution k a = — ^[sign(<j) + -]. The second 
case is £ = — 1 and £' = 1. The solution for this case 
is k b = — ^[sign(q) — -]. Upon a change of variables 
in the energy delta function of Eq. (28) the sum can be 
evaluated. The transfer factor now becomes 



t(Q) = 



r p (q) (M{qj)\ fN (qj)L 



dN 
dT 



\ k B T 2 J \ Airhv 

x [f Q (hvk a ){i - f (hvk b )} 
- f (-hvk b ){i - M-hvk a )}} 



(30) 



where L is the tube length. Substituting the transfer 
factor into our expression for the drag, Eq. (22), gives 



Sd 



\e\Tv 



rag 



aT 2 2nk B hv 



c 3^g n (<l)M{qj)T p (q)N (qj) 



x [f Q (hvk a ){l - f (hvk b )} 

- h{-hvk b ){\ - f (-hvk a )}} 



(31) 



To evaluate the sum over j it is convenient to assume 
that only one linear phonon branch contributes to the 
thermopowcr. The inclusion of other linear modes with 
approximately the same phonon speed will simply add to 
the overall constant. Passing to the continuum limit and 
imposing a Debye cutoff, qo, one can show that 



\eWcL 2 



dqM(q)T p (q)N (q) 



rag ~ 2ir 2 aT 2 k B h J 

x [Jo(^ a ){l-/ (^fc b )} 
-/o(-^){l-/ (-^fc a )}] 



(32) 



To further simplify the drag formula we make a change 
of variables with the following definitions 



hcq (v 



2k B T \c 
1 



e {x-n/k B T) _ l 
V + c 



7 = 



The drag can then be written as 



(33) 



'S'drag — 



\eWcL 2 



rag " TT 2 ah 2 Tv 



2CT du- MT * 



2uc/v ^ 



x [/(-« 7 ){1 - /("«)} - /(«){! - /(«7)}] (34) 

where the Debye temperature is defined as usual Td = 
Hcqn/ks. In the above we have used the fact that for 
(10,10) arm chair carbon nanotubes c = 20.35 x 10 3 
m/s [23] for the longitudinal acoustic mode whereas 
v = 8.4 x 10 5 m/s. This gives f - 100. Equation 
(34) is our primary result and contains several inter- 
esting features. First note that when /i = we have 
f(—x) = 1 — f(x), in which case the two terms in the in- 
tegral cancel giving no contribution to the thermopower. 
If we take \i > then the second term dominates, giving 
a negative contribution to the thermopower. Similarly, 
if we take \x < we get a positive contribution to the 
thermopower. 

The overall sign of Sdrag can be deduced from Fig. 1. 
For fi > 0, intcrband transitions above the crossing 
point are favored. Transitions with q > yield a pos- 
itive change in the electron group velocity which, from 
Eq. (16), give an overall negative sign to the drag. For 
H < 0, interband transitions below the crossing point are 
favored which give and overall positive sign to the drag. 

The Debye temperature plays a small role because the 
kernel of the integral is sharply peaked. As long as the 
peak, u p k — -j^hf , lies within the range of integration the 
Debye temperature can be set to infinity. When the peak 
crosses the range of integration, /i ~ kB ^ v ; 5d rag falls 
to zero. For, [i > kB ^ v the set of transitions favored by 
the electron distribution require wave vectors above the 
Debye cutoff, qo- 

Fig. 2 shows the temperature dependence of 5d rag - To 
evaluate the integral we have taken the deformable ion 
model, Eq. (27). We have also taken a temperature inde- 
pendent phonon relaxation time , m = in Eq. (24). For 
large temperatures Sdrag flattens because we assume here 
that t p does not depend on temperature. In real samples 
it is likely that, for large temperatures, phonon-phonon 



5 



scattering may impose a temperature dependence to the 
phonon relaxation time. We would then see a change 
in the drag at large temperatures. In Fig. 3 we com- 
pare the temperature dependence of the drag for the two 
cases m = and m = 1 in Eq. (24). Here we have taken 
4> = 100K, v/c= 100, n = 0.436 eV, and T D = lOOOK". 
A suppression of the drag at large temperatures clearly 
induces a peak, similar to those observed in experiments. 
The low temperature behavior can be extracted with a 
few approximations. 



Fitting Eq. (2) to measurements of Stotai versus tem- 
perature yields a rough estimate of fi. Fig. 4 shows the 
total thermopower for three possible values of fj,. Mo- 
tivated by experiments on mats of single-walled carbon 
nanotubes [3] we assume that the parameter A changes 
sign with the Fermi level. We have kept the ratio |^| 
fixed. In experiments we expect the parameters A and B 
to have a nontrivial dependence on /i due to sample de- 
pendent variations in the phonon and electronic density 
of states. 



V. THE LOW TEMPERATURE LIMIT 

A low temperature approximation to the integral for- 
mula for Sdrag, Eq. (34), can be obtained for nonzero 
values of fi. For \i < the second term in the integrand 
vanishes. We then have 

\e\rcL 2 
13 dras " n 2 ah 2 Tv 

x / d ^^W- x [7(-«7){l - 7(-«)}] (35) 

As mentioned earlier the factor /(l — /) in the integral is a 
sharply peaked function centered at the value u p k ~ jr^f ■ 
Thus 

_ \e\rcL 2 M(2fj,/hv)T p {2fi/Hv) 
drag ~ ir 2 ah 2 Tv e 2u p" c / v - 1 

x / 2 " d«[/(«){l-7(«7)}] (36) 
Jo 

The integral can now be performed giving 

The above argument remains the same for the case /i > 
except for an overall sign change. The low temperature 
drag contribution to the thermopower becomes the sec- 
ond term in Eq. (2), where the temperature independent 

I I -L 2 

factor B = '^2 a c h2v M(2^/hv)T (2^/hv)(f> m may be taken 
as a fitting parameter. Note that the low temperature 
limit is easily reached because - 3> 1. Fig. 2 compares 
the above approximation to the full formula given by 
Eq. (34). 

Fermi statistics restricts the interband transitions to lie 
near the Fermi level, thereby excluding all but a narrow 
range of phonon wave vectors. Eq. (2) therefore applies 
to systems with any, non-singular q dependence in the 
electron-phonon matrix elements. Only the constant B 
will change with different forms of electron-phonon cou- 
pling. 



VI. CONCLUSION 

We have shown that within a one dimensional model 
phonon drag resulting from interband transitions be- 
tween two linear bands gives a nonlinear temperature de- 
pendence to the thermopower when the Fermi level does 
not lie at the band crossing. Assuming that mechanisms 
other than the electron-phonon interaction contribute to 
phonon scattering we derive a simple expression for the 
phonon drag contribution to the thermoelectric power in 
a model that approximates parameters found in single 
walled (10,10) carbon nanotubes. The strength of the ef- 
fect depends strongly on the position of the Fermi level. 

The above results suggest that the phonon drag effect 
is a good candidate for recent, low temperature anoma- 
lies in thermopower measurements on single walled car- 
bon nanotubes. We note, however, that, as is typical in 
theories of the phonon drag effect in conventional metals, 
the size of the fitting parameter B is difficult to estimate 
from first principles. 

In applying the above formalism to low temperature 
peaks in thermopower measurements under different am- 
bient conditions it is important to account for different 
scattering mechanisms. Different scattering mechanisms 
in nanotubes can exhibit drastically different behavior, 
in analogy to the wide variety of phonon drag effects 
seen in conventional metals with different alloys. In this 
work we have assumed that the electron relaxation time, 
t, is independent of temperature and that t p goes like 
1/T m . From Fig. 4 we see that phonon-phonon scat- 
tering induces a peak in Stotai- The m = 1 case there- 
fore appears be a good approximation for the samples of 
Ref. [3] . Fig. 3 also demonstrates that a temperature in- 
dependent scattering mechanism (m = 0) produces only 
a knee in Stotai- Weak inflections in the thermopower 
have also been observed in many thermopower measure- 
ments on carbon nanotubes. Other conditions may favor 
a decay mechanism which can significantly alter the tem- 
perature dependence of the drag thermopower. Detailed 
comparison with experiment will require suitable choices 
for the temperature dependence of the phonon lifetime. 
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FIG. 1. The linear band model depicting energy along the 
y-axis and wave vector along the x-axis. The two bands 
shown, I = +1 and —1, cross at zero energy. The arrows 
indicate possible interband transitions from an initial elec- 
tron state k to a final state k' . The dashed line indicates 
the position of the Fermi energy relative to the band crossing 
point, defined to be pi. 
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FIG. 2. The phonon drag contribution to the thermopower 
plotted, in units of the constant B, as a function of temper- 
ature. The Fermi energy is chosen to be p = — 0.43eV. The 
ratio between the electron and phonon speeds is taken to be 
- = 100. The solid line shows the result from the full for- 
mula, Eq. (34) with T D = lOOCFT . The dashed line shows 
the simplified approximation, the second term in Eq. (2). For 
purposes of comparison, a temperature independent phonon 
relaxation time is assumed, m = in Eq. (24). 
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FIG. 4. The total thermopower, Eq. (2), plotted versus 
temperature for two values of the Fermi level, fi = — 0.5eV 
and 0.5eV. The ratio of the two fitting parameters is taken 
to be |^4/B| = 10 5 . The top curve has a positive value for 
A while the bottom curve has a negative value for A. We 
also have - = 50, <j> = 100 A', and m = 1. The central line 
is the /i = case in Eq. (34), where A = 0. For /i = 
transitions above and below the Fermi level cancel to give no 
net thermopower. 



FIG. 3. The phonon drag contribution to the thermopower 
plotted versus temperature for two different cases. The 
solid line shows the drag when a temperature indepen- 
dent phonon scattering mechanism dominates phonon decay, 
m — in Eq. (24). The dashed line shows the drag when 
phonon-phonon scattering dominates, m = 1. The Fermi level 
is chosen to be fi — —0.5 eV. The parameter B is the same 
for both curves. We also take ^ = 100 and <f> = WOK. 
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